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Abstract 

We give overcrowding estimates for the Sine/j process, the bulk point process limit 
of the Gaussian /3-ensemble. We show that the probability of having at least n points 
in a fixed interval is given by n —)• oo. We also identify the next 

order term in the exponent if the size of the interval goes to zero. 


1 Introduction 

The Gaussian /3-ensemble is a family of point processes on the real line given by the joint 
density function 

m 

Pm,/3(Ai,...,A,„) = -— Yl |Ai - Ajl^ X (i) 

^m,3 -I ^ ^ -1 

Here /3 > 0 and is an explicitly computable normalizing constant. For the specihc 
values of /3 = 1, 2 and 4 one obtains the joint eigenvalue density of the classical random 
matrix ensembles GOE, GUE and GSE. 

We study the Sine /3 process, the bulk scaling limit of the Gaussian /3-ensemble. Fix 
/3 > 0, and consider a sequence of hnite point processes Am ,/3 with distribution given by the 
density (1). Then as m —)■ cxo the scaled processes 2y/mKm,i3 converge in distribution to a 
translation invariant point process which we call Sine/ 3 . The existence and description of the 
limiting process in the /3 = 1, 2 and 4 cases were obtained by Dyson, Gaudin and Mehta (see 
e.g. the monographs [1], [9], [13]), while in the general /3 case this was done in [14]. We note 
that in our normalization the Sine /3 process has particle density 
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We study the problem of overcrowding for the Sine^? process. We describe the asymptotic 
probability of having at least n points in a hxed interval for a large n. Since the process is 
translation invariant we may assume that the interval is of the form [0, A] with A > 0. We 
will use the notation Np{\) for the counting function of the process, for A > 0 this gives the 
number of points in [0, A]. Our main result is as follows. 

Theorem 1. Fix Aq > 0. Then there is a eonstant c > 0 depending only on {3 and Aq so 
that for any n > 1 and 0 < \ < Xq we have 

g-§n 2 log(^)-cnlog(n+l)log(^)-cn 2 ^ > u) < log(f)+cn 2 ^ ^ 2 ) 

The theorem implies that for a hxed A > 0 we have P(iVa(A) > n) = g^g 

n —)■ oo. If we let n —)> oo and A —)■ 0 simultaneously then we get 

P{Ni 3 {X) > n) = e-|^^^°gO/0+c’0^+»^iogniog(A-i))_ 

Dyson used Coulomb gas arguments in [ 6 ] to predict the asymptotic size of P{Ny{X) = n) 
for a hxed n and large A. His predictions were made precise in [5] for = 1 , 2 and 4 using 
the fact that in these classical cases the probabilities in question have Fredholm determinant 
representations. For general (3 the special case of n = 0 (i.e. the large gap probability) was 
treated in [15]. 

Fogler and Shklovskii in [ 8 ] extended Dyson’s method to give predictions on the asymp¬ 
totic size of P(iVa(A) = n) in the regime where n and A are large. Although this does not 
cover the case we are interested in (A bounded, n large), our results agree with their predic¬ 
tion in the case when n ^ X. Section 14.6 of [9] gives a nice overview of the Coloumb gas 
method in the study of the distribution of A]g(A). 

Our proof uses techniques that we developed in [10], where we proved a large deviation 
principle for ^Ng^X) with scale A^ and a good rate function of the form {31 (p). One of the 
consequences of this large deviation result is that if p > ^ then P(A]a(A) > pX) decays 
asymptotically as as A ^ oo. It was shown in [10] that J(p) grows like 

ip^logp as p ^ cxo so the decay of P(A^^(A) > pA) is formally consistent with our Theorem 
1, even though it describes a different regime of n, A. 

Rigorous overcrowding estimates have been obtained for certain non-Hermitian random 
matrix models. Here one is interested in the asymptotic probability of hnding at least n 
points in a hxed disk. In [11] Krishnapour identihed the overcrowding asymptotics for the 
complex Ginibre ensemble and also gave similar estimates for zeros of certain Gaussian 
analytic functions. In [3] and [2] these results were generalized to point processes obtained 
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by considering products of complex and quaternion Ginibre ensembles. It is interesting to 
note that in all of these non-Hermitian cases the probability of finding at least n points in 
a hxed disk decays like where c is a parameter depending only on the model 

and not the size of the disk. This is of the same form as our result for the Sine^ process. 

The finite version of our problem would be to study the probability of having a large 
number of points in a small interval of order in the bulk of the support of (1). This 

problem has been considered for the spectrum of general Wigner matrices in the context of 
showing the universality of the bulk limit process, see e.g. Theorem 3.5 of [7] or Corollary B.2 
of [4]. These results concentrate on giving sufficiently strong (but not necessarily optimal) 
upper bounds for the overcrowding probabilities. 

The starting point of our proof is the following characterization of the counting function 
of Sine^, due to [14]. 

Theorem 2 ([14]). Let A > 0 and consider the strong solution of the following diffusion: 

dax = X^e~^^dt + 2 sm{ax/2)dBt, aA(0) = 0. (3) 

Then limi^oo exists almost surely, and it has the same distribution as Ng{X). 

This describes the one-dimensional marginals of the counting function Ng{X). One can 
also describe the distribution of the whole function iYg(-) in terms of a one-parameter family 
of diffusions (see [14]), but we will not need this here. The following proposition summarizes 
some of the basic properties of the diffusion (3). Most of these properties are straightforward, 
see Proposition 9 in [14] for some of the proofs. 

Proposition 3. Let A > 0. 

(i) The process is nondecreasing in t. In particular for any integer n we have 

ax{t) > 27rn ift is larger than the first hitting time of27rn. 

(a) //0 < A < A' and ax, ay are strong solutions of (3) with the same Brownian motion 
Bt then ax{t) < ay{t) for all t > 0. 

(Hi) For any integer n the process ax{t) + 2mi solves the same SDE as (3) hut with an initial 
condition of 2 n 7 r instead of 0. 

(iv) For any T > 0 the process fift) = ax{T + t) solves the same SDE as (3) with Ae“ 4 ^ in 
place of X, Bt+T in place of Bt and initial condition 0(0) = ax(T). 
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An important consequence of the previous proposition is the following observation. 

Lemma 4. Fix n > 1 and A > 0. Let t\ he the first hitting time of 271 of the diffusion a^. 
(Note that t\ might he oo). Then 

PiNgiX) >n)=E [^((rA)], where g{t) = P {Nfj{\e~^^) > n - 1 j , 51 ( 00 ) = 0. (4) 

Proof. According to Proposition 3, if we condition on {ta = T} with a T > 0 then the 
process {ax(t + T) — 271,t > 0} will have the same distribution as {a _i3.^{f),t > 0}. Thus 

P[Wm. a\{f) >27m\Tx = T) = P{\im. a _/? (t) > 27r(?7, — 1)). (5) 

The statement of the lemma now follows by taking expectations in tx and using the descrip¬ 
tion of A(g(A) from Theorem 2. □ 

We record the following simple corollary of Lemma 4. 

Proposition 5. For any A > 0 and n > 1 we have 

P(Nf,(X) >n)< , ( 6 ) 

Proof. From (3) one gets Eax{t) = A(1 — e“ 4 *) and using the properties of the diffusion ax 
one can show Elimt^ao ctx{t) = A. Markov’s inequality now implies the statement of the 
proposition for n = 1. For larger values we can use Lemma 4 and the monotonicity of g{t) 
from (4) to get 

P{Ng{X) >n)< P{tx < oo)P{Ng{X) >n-l)< P{Ng{X) > l)P{Ng{X) >n-l) 
from which the statement follows by induction. □ 

Our proof of Theorem 1 will rely on the recursion equation (4) and estimates on P{tx < t). 
We will show that if we can get a good enough upper or lower bound on P(iV^(A) > rr — 1) 
for all A < Ao, then we can also get good bounds for P{Ng{X) > n) for all A < Aq. The next 
section will contain the necessary hitting time estimates while the remaining sections will 
contain the proofs of the upper and lower bounds in Theorem 1. 
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2 Hitting time estimates 


We will estimate the hitting time t\ by coupling the diffusion a\{t) of (3) to a similar diffusion 
with a constant drift term. Let a\ be the strong solution of the SDE 


dax = Xdt + 2 sin ( ^ ) dBt^ '5a(0) = 0 , 


(7) 


using the same Brownian motion as (3). Denote by fx the first hitting time of 27r by ax- 


\" = 


P{tx <t) < P{fx' < t), 

P{rx <t)> P{fx" < t), 

Our bounds will use certain special functions which we now briefly go over. For a < 1 we set 


< ciA'(^) all f > 0 with A' = 

|A, and 

This yields 


with A' = —A, and 

4 

(8) 

with A" = 

4 

(9) 


K(a) = 


t/2 


dx 


t/2 


\/r 


E{a) = / V i — osin^ xdx. 


( 10 ) 


'0 Vi — asm X Jo 

These functions are the complete elliptic integrals of the first and second kind. We note the 
identity 

1 


1 r 


-dt, 


a > —1. 


( 11 ) 


'-00 x/ cosh ^z + a 

We denote by W the lower branch of the Lambert W function (or product log function). 
This is defined as the unique solution of 

^ with W(z) < -1. (12) 

Note that W has domain [—1/e, 0), it is strictly decreasing and satisfies W(xloga;) = logx 
for 0 < a; < 1/e. 

The following proposition summarizes some of the relevant properties of the introduced 
special functions. 

Proposition 6. There exists a constant c> 0, so that the following bounds hold. 

1 


K{-a) 


2 \/a 


log(16a) 


-^^°s(a), \E{-a) - x/d\<-^\og{a), if a > 2. (13) 


\K ^{x)+x ^W^(—a;/4)| < c, if 0 < x < 1/2. (14) 

|>V(-x) - >V(-|/)| < c|log(x/|/)|, |>V(-x)| < clog (i) , if X,//e (0, ^]. (15) 
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Outline of the proof. The estimates (13) are from the Appendix of [10]. The dehnition of 
W allows us to check that the inverse of the function 2 ^ 1 og( 16 a) on [e^/16,oo) is exactly 
—x/4). From this fact and the hrst bound in (13) the proof of (14) is a simple 


exercise. Finally, (15) follows from the observation that W (X) = i □ 


□ 


Now we are ready to state the estimates on the hitting times t\ and fx. 


Proposition 7. There exist positive constants e and c (depending only on ft) so that if 
Xt < e then 



For any Aq > 0 there exist positive constants e,c (depending only on \o,ft) so that if 0 < 


X < Xo, Xt < e and t < elog(^) then 



(17) 


Moreover, the same statements hold for the hitting time fx (with possibly different constants 
e and c). 

We will prove the upper and lower bounds separately. In both cases we will prove the 
estimate for fx hrst, then use ( 8 ) and (9) to get the corresponding statement for tx- 

2.1 Proof of the lower bounds in Proposition 7 

By introducing the new variable Xx = log(tan(Q!A/4:)) we can transform (7) into the SDE 



(18) 


Under this transformation fx becomes the hitting time of oo for the process Xx- Following 
[10] we compare Xx with another family of diffusions. For a > —1 consider the SDE 




and denote by ry,A,a the hitting time of oo. In [10] the Cameron-Martin-Girsanov formula 
was used for diffusions with explosions to compare Xx and Y^.a and the following results 
were obtained. 

Proposition 8 ([10]). For a > —1 and si < S 2 we have 


P{fx e [si,S2]) = E l{TY,X,a ^ [si,S2])e 


( 20 ) 
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Here Gry^aO^) Girsanov factor for which the following estimate holds: 


^TY.X.aO^) + 


X^aTy, 


A,a 


-\{{l + a)K{-a)-E{-a)) 


< 




TY,X,b 


( 21 ) 


The proof of these statements are contained in Section 3 of [10]. Note the slight notational 
inconsistency, in the statements in [10] —a is used in place of a. 

Our lower bound proof will rely on the following lemma, which will be proved at the end 
of this subsection. 


Lemma 9. There is a constant c > 0 so that for any a, X with Xy/a > 2 and a > 2 we have 


P ( Xryxa £ 


«■(-“)(!-a75).«'(-“)(1 + 1^)J) (22) 

Now we are ready to prove the lower bounds in Proposition 7. 

Proof of the lower bounds in Proposition 7. We use (20) and (21) with an a > 0 to get 

P{Xfx < S 2 ) > PiXfx e [si,S 2 ]) = E l(Ary,A,a e 


> E 


l(Ary,A,a e 


Assuming Xy/a > 2 we see that the exponential term under the last expectation is nonde¬ 
creasing in Ty^x^a and because of the indicator function we get the bound 


P{Xfx < S2) > P{XTy,x,a G [5^, S2] . 


Now choose [si,S 2 ] = 4iP(—a)(l — i^), 4iP(—a)(l-|-^^) , and assume a > 2. With 

Lemma 9 we can estimate the probability on the right and we can use the asymptotics (13) 
of Proposition 6 to estimate the exponential term. Collecting all the terms gives 


P(AfA <4A:(-a)(l + 35 ^)) >exp -^K{-a) 


G{1 + A) log(a) 


(23) 


with an absolute constant G. 

Now assume that Aq > 0 is hxed, choose a small e > 0 and assume that 0 < A < Aq, 
AT < i and T < elog(^). Set a = If i is chosen small enough then using the 

bounds of Proposition 6 one can readily show that 


1 iQg(i^) 

2 AT 


<x/d<2 


log(]^) 

AT 


(24) 
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which yields both \^/a > 2 and a > 2 (using A < Aq) for small enough d. Thus we can use 
(23) to get 

logp (ta < r(i + iog|i°[>T) )) + 

>-|^h-f)-ciog(y) 
>-|wy-Ar)-Clog(y) ( 25 ) 


where we again used Proposition 6 and the constant C can change from line to line (but only 
depends on Aq and e). Now set t = T , we can see that T < t < T(1 + lOe). 

The bound (25) gives 


hgP{fx <t)> + 


lOT 


t 


iog(i/(Ar)). 


C'l 0 g(T^ 


2 2 nT 1 

= - f^>v 2 (_AT)-=- 

t ^ ^ ^ ^ ^log(l/(AT)) 


-Cfog(3^) 


Using the second bound of (15) (which is allowed since AT < i) together with 1 < |^ < 1 + lOe 
we get that 


20T 


>V2(-AT)- 


t ' ^log(l/(AT)) ' 

Using the first and then the second bound of (15) we get 


+ Ul 0 g( 3 ;L)<C'fog(fo) 


J-W^-XT) > + Clog If I ■ W(-\t)) > -hv\-\t) - y log(i). 


Putting everything together gives 

fog^^(rA <t)> - c 1^1 + log(ii) 

with a constant c. The only thing left to show is that for small enough e' > 0 if At < e', 
t < £'log(l/(Af)) then we can find T with t = T(1 + iog(i/^r)) )which satisfies AT < e and 
T < elog(^). This is fairly straightforward, one can check that if e' is small enough then 
we can find a T with (1 — <T <t satisfying all the listed conditions. 

Finally, to prove (17) for tx we use the comparison (9). Note that if we have the bound 
A < Ao then A" = |e“ 4 *A < Aq = f Aq. By choosing a small enough e in the bounds \t < e 
and t < e log(^) we can use the just proved lower bound for ta" to give 

fogT(rA <t)> - c(^ + (^fog(if) + f ^ • 

The fact that the right hand side is bounded from below by —|W^(—At) — jlog(^) with 
a new constant c' follows along the line of the already seen arguments using Proposition 6 
and our assumptions on A,t. □ 










Now we return to the proof of Lemma 9. 


Proof of Lemma 9. Our main tool is again a coupling argument. We introduce the diffusions 


A 


dhl = 2 V + a dt — -dt + dB 
dYj^ = — \J cosh^ F+ + a dt + -dt + dB 

driven by the same Brownian motion as Y. Comparing the drifts of Y_, Y^ to that of Y we 
see that Y_ <Y and Y <Y^ until the blowup of the larger diffusions. 

We also introduce Z_ = Y_+t/2 — Bt and Z_|_ = Y^—t/2 — Bt which satisfy the differential 
equations 

Z'_ = —\J cosh^(Z_ — t/2 + Bt) + a > 0, and cosh^(Z_|_ + t/2 + Bt) + a > 0 

until their blowup times. Note that Z_ blows up exactly when Y_ blows up, and Z_|_ blows 
up when Y^ does. Thus Z_|_ will blow up before Y does and Z_ will blow up after. 

We introduce the constants 


£i = 


lOK(-a) 


£2 = 


K(-a), 


6 = 


X^/a ’ ^ 8A 

and we record for further use that since Xy/a > 2 we have 

4^2 + 0 

El > 


—aj T El 
A ^ 


(26) 

We will show that on the event H = {supj<g \Bt\ < E 2 } the process Z_ will blow up before 9 
and the process Z+ will blow up after A“^(4iC(—a) — 2 ei). The final bound will come from 
an estimate on the probability of the event H. 

Note that Z'_ > > 0 until its blowup, in particular Z_ is strictly increasing. Moreover, 

on the event H ioi t < 6 we have 

Z_{t) — 9/2 — E 2 < Z_(t) — t/2 + Bt < Z_{t) + E2- 

Let 7]i be the hitting time of —E 2 by Z_. Then for t < 7]i /\ 9 on the event H we have 

Z'_ (t) = — \J cosh^(Z_(f) — t/2 + Bt) + a > — cosh^(Z_(t) + E 2 ) + a. 

Rearranging, integrating and using identity (11): 

z'_{t) 1 


\ prjiAd 

^iViX9)< 


-.dt < 


cosh^(Z_(t) + E 2 ) + a 


cosh^(z + E 2 ) -L a 


zdz = K{—a), 
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which gives 


\{rii A9) < 2K{—a) < \9, and Xrji < 2K{—a). 


Next we bound the travel time of from —£2 to 6*/2 + £2 which we denote by 772- Since 
ZL > |\/a we have the bound 



Note that we have rji + r]2 < X~^{2K{—a) + £1) < 9 on H. Consider now the travel time 773 
the from 9/2 + 82 to 00 . Denote by 773 the minimum of 773 and 6* — 771 — 772 > 0 . On the event 
H if Z_ > 9/2 + 82, rii + ri2 <t < 9 then we have the lower bound 



Rearranging the inequality and integrating from 771 + 772 to 771 + 772 + 773 we get 

A 773 < 2K{—a). 

Collecting our estimates we see that on the event H we have 

X{rji + 772 + 773) < AK{-a) + Si = X9. 

This immediately implies 773 = 773 and since 771 + 772 + 773 is the blowup time of IT we also get 
that TYXa < 9 on H. 

To show that the blow up time of is greater than A“^( 4 R'(—a) — 28 i) on the event H 
we carry out a similar calculation. The process Z_|_ is strictly increasing and on the event H 
ioi t < 9 we have 


Z^(t) — £2 ^ + t/2 + < Z_{t) + 9/2 + 82 - 


Let Cl be the hitting time of —9/2 — 82 for process Z^. Since the blowup time of Z^ is 
smaller than that of we get that Ci < on the event H. For < —9/2 — 82 on H we 

get 

Z'+{t) = -y^cosh^(Z+ + t/2 + Bt) + a < -^cosh^(Z+(f) - 82 ) + a, 

which gives us 



dt = K{—a) — 




9/2-262 ../cosh^( 2 ;) + a 


1 


dz 


9 + 4:80 

>K{-a)-^^>K{-a)- 8 ,/ 2 . 
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In the last step we used (26). Now take (3 to be the time spent by traveling from £2 to 
00 . Again note that the blow up time of Z_|_ is bounded by that of Z_ which is bounded by 
6 on H. Therefore on H for 82 < Z^ we have 


Z'_^_(t') — — cosh^(.Z_|_(f) + t/2 + Bt) + a < —cosh^(.Z_|_(f) + 6/2 + 82 ) + n. 

This yields 


A 


Cs > K{-a) - 


f‘el2+2e2 


\J cosh^( 2 ;) + a 


dz > K{—a) — 81 / 2 . 


These estimates give 

A(Ci + C3) >4Ja-a)-2£i. 

Our last bound shows that on H the process Y blows up after A“^(4A'(—a) — 28i). Putting 
everything together yields 


P{Xt 


Y,X,a 




-“)(!+3^) )>pm 


On the other hand we have 


P{H) = P{ sup \Bt\ < 82 ) = P(sup \Bt\ < 

t<A-i0 t<l 


and 




A-(-a) 






It is known that P(supi<i \Bt\ < S) > 2(2<h((5/4) — 1) for 5 > 0 (see e.g. Proposition 2.8.10 in 
[12]) which can be bounded from below by cS if we have an upper bound on 6. Since a > 2 
we have an upper bound on ■^x/K{—a) which hnishes the proof. □ 


2.2 Proof of the upper bounds in Proposition 7 

The proof of the upper bound will rely on the following estimate on the moment generating 
function of fx. 

Proposition 10 ([10]). For any a > —1 we have 

^ ^-\{{l+a)K{-a)-E{-a)) ^ ^- 2 ^^ 

This is Proposition 8 in [10], note that in that statement —a is used in place of a. 
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Proof of the upper bound in Proposition 7. By Markov’s inequality we have 


Pifx <t)< {Ee^^^)e-^^ for A<0. 

We will choose A = with a = —K^^^tX/A) > —KP^\e/A) > 0. Using 

Proposition 10 we get 

^-\({l+a)K{-a)-E(-a)) 


This gives the upper bound 

p(^~^ < t) < (,-A{^+AK{-a)-E(-a))+P^t+t^t ^ ^K^-i){t\/i)^+\E{-a)+^ ^ 


where we used our specihc choice of a. If e is small enough we will have a > —KP^\e/A) > 2 
which allows us to use the asymptotics of Lemma 6. This gives the bound a < log^(^) 

with a positive constant U, and also 


KP^\tX/A)XH 

8 


+ AU(— cl) + 



< -lyvi-XtHf + ^ + C(1 + i) log(i) 


(with a possibly different C). Using the bounds in (15) and changing the value of C again 
we get 


KP^\tX/A)XH 

8 


+ \E{-a) + 



<_|w(_A()^ + C(l + l)log(i), 


which finishes the proof for fx. 

To get the similar statement for tx we use (8) to obtain 

p{tx <t)< P{fi^ <t)< (“T^) ^ ^ (^) ■ 

Note that we are allowed to use the just proved bound on f if we choose the upper bound e 
on At small enough. To complete the proof we use (15) to replace ^ with A everywhere by 
choosing e small enough and modifying the constant in front of the last term of (28). □ 


3 Proof of the lower bound in Theorem 1 

Our proof will have two main steps. In the first step we will give a reasonably good (but not 
optimal) lower bound on P(A(a(A) > n). In the second step we will show how one can get 
better and better bounds as we increase n. 
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Proof of the lower hound in Theorem 1. We first give a direct estimate on P{Nj^{\) > n) for 
A < Aq. Although this lower bound will be far from the statement of the theorem, we will 
later show how to ‘bootstrap’ this into a near optimal bound for large n. 

Recall Theorem 2 and Proposition 3. The connection between the counting function 
Np[\) and the diffusion a\it) shows that for any T > 0 we have 

P{Np{X) >n)> P{ax{T) > 2mr). 

The drift in (3) on [0,T] is at least A' = |e“^^A so coupling ax and ay (using the same 
driving noise) we get 

P{ax{T) > 2n7r) > P{ay{T) > 2n7r). 

For the process ay denote the hitting time of 27ik by fy,k (note that fy,i = fy using our 
previous notation). Using the same ideas needed for Proposition 3 it is easy to show that 
ta',!, fx \2 — ■ ■ ■, ^\',k — ^\',k-i, ■ ■ ■ are i.i.d. random variables. Using this we can obtain 

the estimate 


y[ay{l) > 


= P{ry < 


T\n 


We will set T = ^log(^) and apply the lower bound (17) of Proposition 7 with t ^ and 
A'. (Note that A' < f Aq.) We first have to check that the conditions of the proposition are 
satisfied. From the definitions we get 

t 1 41og(f) 

log(^) n 2/31og(f)-/31oglog(f)' 


= (^) log(f). 


Since we have A' < |Ao, there is a positive integer no (depending only on Aq, /?) so that if 
n > no then both of these expressions are bounded by e. (Here £ is the constant from the 
statement of Proposition 7 which also only depends on Aq,/?.) Thus if n is large enough we 
can use Proposition 7 to get 


P{fy <l)> exp 


- c + ^ ) log(3A) 


= exp 


(dn 


[ 21og(f) 




^iog(g; 


c 1 + 


jSn ^ f n'^ 1 


Choosing n large enough we can use the bounds (15) of Proposition 6 to get 


P{r\' exp [-cinlog(f) - C2n] 
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with positive constants Ci,C2 (which only depend on Pntting onr estimates together 

we obtain the lower bonnd 


P{Nf}{X) >n)> exp log(f) - , (29) 

which holds for n > Uq with uq only depending on Ao,/5. Finally, since P{Np{\) > k) is 
decreasing in k, we can adjust the constants ci, ca so that (29) remains valid for all n > 1. 

To finish the proof we will show how to improve this lower bound to get the lower bound 
in Theorem 1 . From the recursion equation of Lemma 4 we get that for any T > 0 and n > 1 


P{Nis{\)>n)>E[l{T,<T)g{T,)] 

> P{rx < T)g{T) = P{tx < T)P(7V^(Ae-^'^) > n - 1), (30) 


where we used the fact that the g{t) function is non increasing. Set T = —Using the 
same arguments as in the first part of this proof we can check that there is a positive integer 
riQ (depending only on Xo,(3) so that for n > no we can use the lower bound of Proposition 
7 with t = T. This gives 


logP(rA < T) > 


f3n 

2bi(f) 





■bog(x) 


n 


c(l + 


f5n 

41 ^ 


log 



1 

log(f) 


Assuming that n is large enough we can use the bounds in (15) to remove the constant | 
from inside W. This allows us to hnd positive constants ci and ca so that 


logP(rA <T)> 


(3n 


2 \ogC- 




-bog(f 


n 


)^ - cilog(f) - Can 




= --nlog(f) - Cl log(f) - Can 


(31) 


where we used >V(a:logx) = logx which holds for 0 < x < 1/e. 

Now assume that for a certain n there is a constant fn-i > 0 so that we have the lower 
bound 


P(7V^(A) > n 


1) > exp 


-fn-l log(^^^^y^) 


c(n — 1)^ 


for all A < Ao with a fixed constant c. Then 


logP(iV^(Ae 4^) > n - 1) > -/„_i log(^^) - c(n - 1)^ 

Ae 

71 — 1 

= -/„_! log(f) - /„_! log(-) - ^T/„_i - c(n - 1)2 

n 4 

^ n + 1 - In / nn 2 

>- fn-i log(-) - c(n - 1)2. 

n n 
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Together with (30) and (31) this yields the lower bound 


logP(iV^(A) >n) > 


n + l , , 

- tn-l + 77^ + Cl 

n 1 


log(f) 


C2n 


c[n 


1 )^- 


If we also assume that c > C2 then we can further simplify this to 

\ogP{Nfi{X) >n)>- + ^ri + ci^ log(f) - cn^. 

Thus we proved the following statement. Suppose that no is large enough, c > C2 and / 
solves the following recursion: 


Jn — - Jn-l + —n + Cl, 

n 1 


for n > uq. 


(32) 


Then if 

logP(iV^(A) > n) > -/nlog(f) - cn^ 

holds for n = no then it also holds for any n > no- Because of (29), we can choose /„g so 
that the previous inequality holds for no (with a suitable c). To hnish the proof, we need to 
show that the recursion (32) will have the appropriate asymptotics. The general solution of 
the recursion (without the initial condition) is 


fn = ^n{n + 1) + (ci - |) ^ + C{n + 1). 

k=l 

From this it is clear that there is a solution with any initial condition, and the solution 
satishes 

fn < + bn\og{n + 1) 

with a hnite b. From this the lower bound of Theorem 1 follows for n > uq and A < Aq. 
Using the argument following (29) we can now modify the constant b so that the lower bound 
holds for all n which completes the proof. □ 


4 Proof of the upper bound in Theorem 1 

The upper bound will use a two step strategy similar to that of the lower bound, but our 
recursive step will be a bit more involved. 

Proof of the upper hound in Theorem 1. Recall from Proposition 5 that 

P{Np{\) >n)< < exp(-nlog(l/A)). 
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By choosing Ci > 1 we get —nlog(l/A) < —nlog(n/A) + Cin^ which yields 

logP(A^^(A) >n)< —nlog(n/A) + Ciu^. (33) 

Now assume that n>2 and there are positive constants fn-i < f — 1)^, Ci so that for 
all 0 < A < Ao we have 

P{Nf 3 {X) > n - 1) < exp [-fn-i log (^) + Ci(n - 1)^] . (34) 

Our goal is to show that we can choose an uq (depending only on Xq, (3) so that if n > no 
and Cl is larger than a hxed constant then we can get an upper bound of a similar form for 
P{Ni 3 {X) > n) as well. 

Let To = 0 < Ti < • • • < Tm- Then 

M-l 

P{Nfj{X) > n) = ^ P{N/ 3 {X) >n,Ti<Tx< T^+i) + P(iV^(A) >n,Tx> Tm). (35) 

i=0 

Choose real numbers 71 < • • • < 7 m with 71 = 1/3,'Jm = 3, M < 3n and 7^+1 — 7^ < 1/n. 
We set Ti = for i = 1,, M, if n > Aq then these are positive numbers. We will 

estimate the sum (35) term by term. 

Using the identity (5) from the proof of Lemma 4 we can show that 

PiNjsiX) >n,Ti<T< Ti+i) < P(iV^(Ae"4^') > n - l)P(r < T^+i). (36) 


From the dehnition of T* (if n > max(l, Aq)) we get 


(1 + 7A) log(n/A) - < log (^) . (37) 

Note that AT* < ^^log(n/A) which can be made as small we wish by choosing no large 
enough. Thus we can apply the upper bound in Proposition 7 to get 


logP(r < Ti+i) < 


< 


< 


Tj+i 


W^(—ATj+i) + c ( 1 + 


T 


i+l 


log( 


XTi 


+ 1 ' 


n 


P' (")) + C (1 + 


p 


n 


27i+i log(n/A) 

■—^n log(f) + C log(f) + Cn. 
27i+i 


47i+i log(n/A) 


log(f) 


(38) 


Here we used the bounds (15), 7*+! > 1/3 and the identity W^xlogx) = logx. The constant 
C depends only on 13 and Ao, but may change from line to line. 
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Using the assumption (34) for Ae 4 ^^ < Aq we get 


P(7V^(Ae-4^0 > n - 1) < exp -f^-i log ) + c^{n - 1)^ 

V Ae T O , 


Using (37) we can write 


-fn-i log ( ^ < -/n-i(l + —) log(n/A) + 

VAe-4^V n n-1 

< - (^/n-i(l + ^) - 0 log(n/A) + - 1), 

where we used fn-i < §{n — 1)^ and 7j_|_i — 7* < which holds for f > 1. We also have 


log 


AP 


i+l 


= log(n/A) - loglog(n/A) - log(47i+i//3) < log(n/A), 


where we used 1/3 < 7* < 3 and that n/X can be made sufficiently large by choosing hq 
appropriately. 

Collecting all our estimates we get that for 1 < i < M — 1 we have 


logP(iVa(A) >n,Ti<T< T^+i) 


< - ( /n-i(l + ^) + ) log(’^/-^) + ci(n - 1)^ + C'log(f) + Cn. 

n 27i+i 


Assuming ci > C this yields 


log P{Nfi{X) >n,Ti<T< Ti+i) 

< - f/n-i(l + ^^) + - c\ \og{n/X) + cin^ - cin 

V n 27i+i J 

< - login/ A) + cin^ - Ciu. ( 39 ) 

In the last step we used that x(l + ^) + > a: + \/2j3x if a: > 0,7 > 0. 

For the i = 0 case we can use a similar argument (and the fact that 71 = 1) to show that 

log P(i\(g(A) > n,To < T <Ti) < -{fn-i + ^^n - C) log(n/A) + cin^ - ciu. 

Since /„_i < ^{n — 1)^ we have |3 ?t, > \/2f3fn-i and so inequality (39) holds for i = 0 as 

well. 

To bound the last term in (35) we start with the bound 

PiN^W >n,TM<T)< P(iV^(Ae-^"’-) > n - 1). 
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Recalling that 7 m = 3 and using arguments similar to what we used to obtain (39) we get 


P{Nf}{X) >n,TM<T)< -fn-i ( 1 + - ) log(n/A) + cin - cin. 


(40) 


Applying the upper bounds (39) and (40) on the respective terms of the sum in (35) yields 

log P{Nfi{X) >n)<- min + A/2/5/n-i - C 2 , fn-i (^1 + ^ log(n/A) 

+ Cin^ — Cin + log( 3 ? 7 , + 1), 

with constants ci, C 2 that only depend on Aq, /9. By choosing ci > 2 we can simplify this as 


log P{Nf}{X) >n)<- min - C 2 , fn-i (^1 + ^ log(n/A) + 


cin^. 


Combining this result with (33) we get the following statement. There is a positive integer 
Uq, and positive constants c, ci (all depending only on Ao,/S) so that if we set fn^ = uq for 
an no > u'q and 


then 


fn = min + \/ 2 / 5 /n-l - C 2 , fn-1 , for u > no, 


log P{Np{X) >n)< -/nlog (y) + Cin^ 


(41) 


for all n > no, X < Aq. 

To hnish our proof we just have to show that we can choose a large enough no as our 
starting point so that the solution of the recursion (41) satishes /„ > — cnlog(n + 1) 

with a constant c. 

This is fairly straightforward although somewhat technical. A simple computation shows 
that if we choose no large enough then the following statements hold: 

(i) fn is strictly increasing, 

(ii) fn = fn-i (l + f) for no < n < ni for a certain ni, 

(hi) fn = min (^fn-i + y/2/3/„_i - C 2 , for n > ni, and 

(iv) if n is large enough and fn-i > |(n — 1 )^ — c(n — 1 ) logn with a large enough c then 


we also have 


fn-i + \/2/3/„_i - c > - cnlog(n + 1 ). 
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Putting these statements together we get that /„ > — cnlog(n + 1) for a large enough 

n which gives 

log P{Ni 3 {\) >n)< — — cnlog{n + 1)^ log (^) + ciu^ (42) 

for a large enough n. Using (33) we can now modify c, Ci so that (42) holds for all n and 
A < Ao, which hnishes the proof. □ 
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